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ABSTRACT 

This paper reports the first study of the O VI resonance hne emission (AA 
1032, 1038) originating in the Local Bubble (or Local Hot Bubble) surrounding 
the solar neighborhood. In spite of the fact that O VI absorption within the Local 
Bubble has been observed, no resonance line emission was detected during our 230 
ksec Far Ultraviolet Spectroscopic Explorer observation of a "shadowing" filament 
in the southern Galactic hemisphere. As a result, tight 2 sigma upper limits are 
set on the intensities in the 1032 and 1038 A emission lines: 500 and 530 photons 
cm~^ s~^ sr~^, respectively. These values place strict constraints on models and 
simulations. They suggest that the O Vl-bearing plasma and the X-ray emissive 
plasma reside in distinct regions of the Local Bubble and are not mixed in a single 
plasma, whether in equilibrium with T ~ 10^ K or highly overionized with T ~ 4 
to 6 X 10^ K. If the line of sight intersects multiple cool clouds within the Local 
Bubble, then the results also suggest that the hot /cool transition zones differ from 
those in current simulations. With these intensity upper limits, we establish limits 
on the electron density, thermal pressure, pathlength, and cooling timescale of 
the O Vl-bearing plasma in the Local Bubble. Furthermore, the intensity of O VI 
resonance line doublet photons originating in the Galactic thick disk and halo 
is determined (3500 to 4300 photons cm^^ s^^ sr^^), and the electron density, 
thermal pressure, pathlength, and cooling timescale of its O Vl-bearing plasma 
are calculated. The pressure in the Galactic halo's O Vl-bearing plasma (3100 to 
3800 K cm~^) agrees with model predictions for the total pressure in the thick 
disk/lower halo. We also report the results of searches for the emission signatures 
of interstellar C I, C II, C III, N I, N II, N III, Mg II, Si II, S II, S III, S IV, S VI, 
Fe II, and Fe III. 



Subject headings: Galaxy: general — Galaxy: halo — ISM: general — ISM: 
individual (Local Bubble) — ultraviolet: ISM 
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1. Introduction 

The Local Bubble (LB) or Local Hot Bubble and its environment are sketched in Fig- 
ure 1. The Local Bubble, a ~ 10^ pc^ region of X-ray emissive, presumably hot (~ 10^ K) 
plasma, is nestled within a rarefied cavity in the Galactic disk called the Local Cavity, 
(Knapp 1975; McCammon et al. 1983; Snowden et al. 1998; Sfeir et al. 1999). Within the 
Local Bubble, lie a number of cool (~ 10^ K) clouds, including the complex of parsec-scale 
clouds surrounding the Sun (Prisch 1986; Lallement & Bertin 1992; Gry et al. 1995). Esti- 
mates of the number of clouds per line of sight range from ~ 2 to ~ 6, with ~ 2 coming 
from the ratio of the typical observed absorbing column density (Hutchinson et al. 1998) to 
the absorbing column density of the local cloud (Lallement et al. 1995), and ~ 6 coming 
from the number of clouds along the e CMa line of sight (Gry et al. 1995). Beyond the Local 
Bubble, in the direction of the Galactic center, lies a superbubble named Loop I, which was 
blown by the stars and supernovae in the Sco Cen association (Egger 1993). 

The Local Bubble is larger and more energetic than a supernova remnant, yet smaller 
and less energetic than a superbubble. Such bubbles are important as most of the hot 
gas in the Galactic disk resides within them. Like the distribution of small versus large 
stellar associations, the population of bubbles is likely to be larger than the population of 
superbubbles. Nonetheless, because its softer X-ray photons from its relatively cooler plasma 
are more easily absorbed, a bubble like the Local Bubble would be much more difficult to 
detect from a great distance than would an energetic superbubble like Loop I. Thus, many 
other Local Bubble analogs may reside within the Galaxy, unknown to us because they are 
obscured by the neutral and molecular material of the Galactic disk. In this sense, we are 
fortunate to have a local specimen to examine. 

At this time, our understanding of the Local Bubble is chiefiy phenomenological. To 
some degree, the Local Bubble is defined as the source of soft X-rays (~ | keV) produced 
between us and the nearest opaque material. Its physical presence has been surmised from 
the anti-correlation between soft X-ray intensities and neutral hydrogen column densities 
(associated with the larger Local Cavity), detections of X-rays originating in the foreground 
of opaque clouds, relative constancies between Be (0.07 to 0.11 keV), B (0.13 to 0.19 keV) 
and C (0.16 to 0.28 keV) band surface brightnesses across the sky (implying very little 
absorption as the effective absorption cross sections vary significantly between the bands), 
and consistency among two decades of X-ray observations. 

Fundamental characteristics of the LB, such as its presumed temperature and size, 
cannot be determined directly, but can only be estimated when additional constraints (such 
as ionizational equilibrium, constancy of emissivity, and filling factor) are assumed. Given 
these assumptions, the estimated temperature has been derived from the ratio of fiuxes 
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in broad soft X-ray bands. In directions where the Galactic disk blocks all of the | keV 
emission from beyond the Local Bubble, the observed ROSAT R2 to Rl band ratio is ~ 
1.06, implying a temperature of 10^'^ K (Kuntz & Snowdcn 2000). This temperature is 
consistent with values found by shadowing observations (described below) at higher Galactic 
latitudes (Kuntz & Snowden 2000; Snowden et al. 2000) and with Wisconsin's Ultrasoft X-ray 
Telescope (UXT) and All-Sky Survey's B/Be and C/B band ratios (Bloch et al. 1986; Juda 
et al. 1991). Similar temperatures were also derived from the Diffuse X-ray Spectrometer's 
soft X-ray spectra, although these high resolution spectra were difficult to fit with simple 
models (Sanders et al. 2001). 

Desirous of a better understanding of our environs as well as a model of the bubble 
population (useful for calculating the hot gas filling factor), both observers and theorists 
have long sought explanations of the LB's origin and life history. Many researchers have 
pointed out that the energy embodied in the Local Bubble's X-ray emitting plasma is so 
great as to require a supernova or a series of supcrnovac, possibly aided by stellar winds. 
The energy may have been deposited within an existing gap in the H I disk or may have 
pushed aside the disk material to create a cavity. Alternatively, the explosion(s) and winds 
may have originated in the nearby Loop 1 (Sco-Cen) superbubble and blown into the solar 
region (Frisch 1981; Bochkarev 1987), or the Local Bubble and Loop I may have been blown 
independently but are now interacting (Egger & Aschenbach 1995). In some of these scenar- 
ios, the X-ray emitting plasma may have expanded after being shock heated, vastly cooling 
the plasma, retarding the coUisional recombination rates, and thus "freezing" the plasma into 
an "overionized" state (Breitschwerdt & Schmutzler 1994; Breitschwerdt 2001). In contrast, 
if the expansion was restrained and the time-scale was long, the plasma may have already 
approached collisional ionizational equilibrium (Smith & Cox 2001). These two cases should 
produce markedly different spectral signatures. Presently, many theorists and observers arc 
working to decipher the Local Bubble's life story by comparing its X-ray spectra to spectra 
predicted by hydrodynamical simulations. It is hoped that when the observed fluxes in mul- 
tiple bands or emission lines are matched by the predictions of a particular type of model, 
then we will have found a key to deciphering the story of the Local Bubble. Here, we extend 
this sort of comparison into the ultraviolet, using the O VI resonance line (1032, 1038 A) 
intensities. We also compare the intensities with O VI absorption line column densities taken 
from the literature, in order to estimate the physical properties of the O Vl-bearing plasma. 

We use the "shadowing" strategy long employed in X-ray analyses to isolate the Local 
Bubble's emission from the more distant interstellar medium's emission (see Section 2). The 
observations (228 ksec toward I = 278.6°, b = —45.3°) and data reduction techniques are 
discussed in Section 3. Although we observed for a very long time, no O VI resonance line 
emission was detected. Tight upper limits were placed on the 1032 and 1038 A lines (see 
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Section 4). The results of our searches for other cosmic hues within the bandpass are also 
reported (see Section 4). Our upper limit on the O VI doublet intensity is well below the 
4 X 10^ K, overionized model (Breitschwerdt 2001) predictions. This discrepancy, as well as 
discrepancies in the C III intensities and the O VI column density, practically disallow models 
of this sort. Our two sigma upper limit on the O VI doublet intensity is only marginally 
above the minimum intensity expected from the O Vl-rich zone around the most local cool 
cloud (Slavin 1989) combined with the O Vl-rich regions in Smith & Cox (2001) 's suite 
of simulations of a multiple supernovae induced hot bubble. In the event that the line of 
sight intersects multiple cool clouds within the Local Bubble, as well as the Local Bubble's 
wall, then our doublet upper limit is well below the predicted intensity and places a strong 
constraint on models of hot gas/cool gas transition zones. In order to determine if the upper 
limits require the Local Bubble to be undergoing unusual physical processes or contain 
unusual plasmas, we constructed and compared with a generalized multi-component (an 
O Vl-rich ~ 3 X 10^ K plasma and a soft X-ray emitting ~ 10® K plasma) model of the Local 
Bubble. This model's O VI intensity and column density and soft X-ray surface brightness 
were all consistent with observations (see Section 5). By combining our O VI intensity upper 
limit with other researcher's O VI column density estimates, we were able to calculate limits 
on the electron density, thermal pressure, pathlength, and cooling timescale (also see Section 
5). By subtracting the Local Bubble's O VI doublet intensity upper limit from measurements 
for longer, high latitude lines of sights, we found the halo's intensity, from which its electron 
density, thermal pressure, pathlength, and cooling rate were estimated. The results of this 
project are summarized in Section 6. 



2. Shadowing Strategy 

The intensity of resonance line photons from O VI along extended lines of sight through 
the Local Bubble, the Galactic halo, and intervening regions is reported elsewhere (Shelton 
et al. 2001; Dixon, et al. 2001; Shehon 2002; Welsh et al 2002). Here, we wish to measure 
only the O VI intensity of the Local Bubble, and so must block the photons from exter- 
nal regions. To do so, we employ the "shadowing" strategy oft-used in X-ray analyses. In 
the ideal shadowing observation, an opaque object blocks the emission originating beyond 
it, so that the observed emission originates between the viewer and the object. For the 
shadowing object, we have chosen a filament in the Southern Galactic hemisphere (see Fig- 
ure 2). Penprase et al. (1998) estimated the distance to the filament as 230±30 pc, which 
places the filament between the Local Bubble and the Galactic halo. Penprase et al. (1998) 
also estimated the mean color excess {E{B — V)) of the filament as 0.17±0.05 magnitudes, 
implying an absorbing column density of Nh = 8.4 ± 2.5 x 10^° cm~^. This is consistent 
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with the column density imphed by the Schlegel, Finkbeiner, & Davis (1998) IRAS 100 //m 
measurement of ~ 7.3 MJy sr~^, which when scaled by the 100 fim to Nh converstion re- 
lation for the southern hemisphere listed in Snowden et al. (2000), yields a column density 
of 9.9 X 10^*^ cm~^. For this color excess and the extinction curve and parameterization 
presented in Fitzpatrick (1999), 89lf]^% of the O VI photons originating beyond the filament 
should be blocked. 

As a demonstration of the filament's "shadowing" function, we note it's ability to block 
soft X-rays. (The theoretical extinction of ~ ^ X-rays is larger than that of ultraviolet 
photons.) The ROSAT ^ kcV band image (sec Figure 2) reveals an obvious soft X-ray 
depression coincident with the filament. The surface brightness seen in the direction of the 
filament and attributed to the Local Bubble is 590 x 10~^ counts s"^ arcmin"^, while the 
"off-filament" surface brightness is ~ 1300 x 10~^ counts arcmin"^ (K. D. Kuntz, private 
communication). Not surprisingly, the filament has already been used successfully in ROSAT 
soft X-ray shadowing analyses (Wang & Yu 1995; Snowden et al. 2000). 

One potential difficulty in using this region of the sky is posed by Penprase et al. 
(1998) 's suggestion that the filament may be part of a supernova remnant. If this were the 
case, then the foreground O VI intensity should be attributed to both the Local Bubble 

and the supposed supernova remnant. Because the results or our study are upper limits, 
this point is unimportant. Observations of the targeted region, the bulbous portion of the 
filament around / = 278.6°, h = —45.3°, are discussed in the following section. 

3. Observations and Data Reduction 

The FUSE focal plane assembly contains four spectrograph entrance apertures. The 
largest, the low resolution (LWRS) aperture measuring 30" x 30", is used for observations 
of diffuse sources. Furthermore, the FUSE instrument includes multiple approximately co- 
aligned channels. Four of the channels (LiF lA, SiC lA, LiF 2B, and SiC 2B) cover the 
1032, 1038 A region. Of these, the LiF lA is most sensitive and least prone to scattered 
light problems in the 1032, 1038 A region. We use this channel for our O VI analysis. In 
addition, we use the LiF lA, SiC IB, LiF 2A, and SiC 2A channels for the C I, C II, C III, 
N I, N II, N III, Mg II, Si II, S II, S III, S IV, S VI, Fe II, and Fe III analyses. 

The raw data consists of multiple observations of three nearly coincident directions 
(see Table 1 and Figure 3). The combined exposure time (after subtracting time when 
the spacecraft was exposed to the Earth's South Atlantic Anomaly, time when the detectors 
recorded anomalous bursts, etc.) for the LiF lA observations is 228 ksec. The other detectors 
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have had differing useful exposure times. Most notably, detector 2 malfunctioned during the 
12050501 observation, reducing the total observation time with the LiF 2A and SiC 2A 
channels. The observations can be grouped according to pointing direction. Henceforth, the 
combined 12050501 and 11050510 observations of / = 278.58°, b = -45.31° will be denoted by 
"120505", the 12050601 observation of / = 278.59°, b = -45.30° will be denoted by "120506", 
and the combined B1290101 and B1290102 observations of / = 278.63°, b = -45.31° wiU 
be denoted by "B12901". Each of these three datasets was processed using version 2.1 
of the CALFUSE pipehne (Sahnow et al. 2000; Dixon, Kruk & Murphy 2002). We set the 
appropriate pipeline switches, so as to 1.) exclude counts that produced pulseheights outside 
the expected range for cosmic photons (the chosen CALFUSE pipeline cutoffs were 4 and 12 
in the standard arbitrary units for the older datasets (120505 and 120506) and 4 to 31 for the 
newer datasct (B12901)) in order to reduce the detector background noise, and 2.) disallow 
automatic background subtraction, enabling us to more accurately calculate the signal to 
noise of the results. 

Figure 4 displays the 1020 to 1045 A region of the LiF lA spectra for the three pointing 
directions. At this stage, the absolute wavelength scales are still inaccurate, though the 
relative wavelength scales of the pipeline spectra are accurate to several detector pixels 
(about 0.035 A in the LiF lA spectrum). The largest sources of systematic error in the 
following intensity calculations are the uncertainties in the solid angle of the aperture and 
effective area calibration. Both the LWRS's sohd angle (30 arcsec x 30 arcsec) and the 
LiF lA effective area cahbration are thought to be accurately known to ~ 10%. Thus, the 
systematic uncertainty in intensity measurements will be ~ 14%. 



4. Spectral Analysis 

In order to determine the absolute wavelength scales for the LiF lA spectra from the 
three pointings, we used the airglow emission lines as reference wavelengths. We then co- 
added the spectra from the three pointings and shifted to the Local Standard of Rest (LSR) 
reference frame. The 1020 to 1045 A region of the resulting sateUite-night and day + night 
spectra are displayed in Figure 5. 

We found the continuum level in the satellite-night spectrum by fitting a smooth curve 
to the surrounding spectral region. We then searched the residual spectrum for emission 

features. None were found in close proximity to the O VI resonance line rest wavelengths in 
the LSR reference frame. If cosmic features had existed, their widths would have equaled or 
exceeded the instrumental width function, their signals would have been calculated from the 
numbers of counts in excess of the continuum within the extraction width, and their random 
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uncertainties would have been calculated from the square roots of the numbers of spectral 
counts (not residual counts) within the extraction width. We have followed this method, 
using the total instrumental width (generously set to 0.43 A) as the extraction width, to 
determine the upper limits. The "signal" and 1 sigma statistical uncertainty intensities for 
the O VI resonance lines at their rest wavelengths in the LSR reference frame are 180 ± 330 
and -180 ± 310 photons cm~^ s^^ sr^\ for the 1032 and 1038 A lines, respectively. These 
measurements are subject to the above mentioned systematic uncertainties of 14%. We 
validated our technique by applying the following alternate extraction technique to sample 
spectra. We fit the spectrum to a continuum and small emission feature, then varied the 
size of the feature to determine the 1 sigma uncertainties. The results were similar to those 
found by applying our techniques to the same spectra. 

When the sateUite-day portions of the data sets were added, the spectrum became 
more ragged. Prom the variations between the 120505, 120506, and B12901 satellite-day 
spectra we surmise that the raggedness is primarily due to scattered light. For the combined 
day and night spectrum, the "signal" and 1 sigma statistical uncertainty intensities are 20 
± 230 and 80 ± 210 photons cm~^ s^^ sr^^, for the 1032 and 1038 A lines, respectively. 
The measurements are subject to the 14% systematic uncertainties. For reference. Table 2 
tabulates these results and lists the tightest 1 and 2 sigma upper limits. The 1 and 2 sigma 
upper limits were calculated from the unrounded numbers. 

We also searched for C I, C II, C III, N I, N II, N III, Mg II, Si II, S II, S III, S IV, 
S VI, Fc II, and Fc III emission lines in the spectra taken with the LiF lA, SiC IB, LiF 2A, 
and SiC 2A detector segments. We relied on the satellite-night portion of the data when 
we searched for emission lines which also appear strongly in the solar spectrum (i.e. C III, 
Curdt et al. (1997)), appear weakly in the terrestrial airglow spectrum (i.e. N II, Feldman 
et al. (2001)), or are significantly stronger in the sateUite-day portion of the data than in 
the satellite-night portion of the data (i.e. C II). The only signals detected were marginally 
significant intensities of C III and N II (see Table 3). 



5. Discussion 

These observations are the first to place tight upper hmits on the intensity of the O VI 
emission from the Local Bubble (2 sigma upper hmits on the 1032 and 1038 A emission lines 
of 530 and 500 photons cm~^ s~^ sr~^, respectively). The theoretical ratio of the intensity 
of 1032 to 1038 A photons is 2.0 to 1.0. Using the 2.0 to 1.0 ratio and the 1032 A line upper 
limit yields a 2 sigma doublet upper limit of 800 photons cm~^ s~^ sr~^ or 1.5 x 10~^ ergs 

—2 —1 —1 

cm s sr . 
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5.1. Compcirisons with Local Bubble Models 

Breakout Model: 4 x 10"^ K, Overionized Bubble 

In this scenario, a series of supernova explosions and winds has created a hot, youthful 
bubble. It is assumed that the bubble expanded rapidly as it broke into a low density region. 
The rapid, adiabatic expansion cooled the bubble faster than the ions could recombine, so 
that the gas now has a temperture of 4 to 6 x 10^ K, but still highly ionized and X-ray 
emissive (Breitschwerdt & Schmutzler 1994; Breitschwerdt 2001). Because the transition 
zones between the bubble plasma and the embedded cool clouds will be cooler than the 
coUisional ionizational equilibrium temperature for O VI, they need not be rich in O VI ions. 

In the Breitschwerdt (2001) model, the Local Bubble has a radius of 115 pc, density 
(of electrons, Breitschwerdt, personal communication) of 2.4 x 10~^ cm~^, and pressure 
(presumably thermal) of ~ 2000 K cm~^. Thus, the plasma temperature is approximately 
4.3 x 10^ K. Lacking published estimates for the O VI intensity, we will calculate it from the 
above values and the emission equation in ShuU & Slavin (1994). The intensity calculation 
also requires the ratio of electrons to hydrogen nuclei in a solar abundance plasma (1.2), the 
solar abundance of oxygen (8.5 x 10~^, Grevesse & Anders (1989)), and an estimate of the 
fraction of oxygen atoms in the O VI ionization level. During the plasma's extraordinary 
cooling history, it has experienced some degree of recombination, as evidenced by the strong 
lines of C III, O III, and O IV in Breitschwerdt (2001) 's predictions. Here, we take the 
lower limit on recombinations of oxygen to the O VI state as that of a 5 x 10^ K coUisional 
ionizational equilibrium plasma. In this case, the fraction of oxygen in the O VI state is 
0.037 (Mazzotta et al. 1998; Schmutzler & Tscharnuter 1993). As the plasma recombines 
further, the O Vl/oxygen fraction increases, and, due to recombinations of O VII, does not 
decrease until the most of the oxygen has recombined to O I and O II. An upper limit on the 
extent of the recombinations is taken from the isochoric cooling predictions of Schmutzler & 
Tscharnuter (1993) at the model temperature, again yielding an O VI to oxygen fraction of 
0.037. Wc take 0.037 as our estimate, noting that the O Vl/oxygen ratio and the intensity 
prediction may be higher. The resulting intensity estimate, ~ 1900 photons cm^^ s^^ sr^^, 
is more than twice the 2 sigma upper limit established in this paper. 

For this calculation, the O VI ions are assumed to be mixed throughout the remnant. 
If the assumed abundance of metal atoms were reduced from the solar value, then the soft 
X-ray scaling relationships used in finding the model parameters would need to be adjusted 
to a similar degree. As a result, the model parameters would need to be adjusted, probably 
leading to a bigger, denser, and/or hotter bubble. As far as the predicted O VI intensity is 
concerned, these adjustments would probably offset the hypothetical abundance adjustment. 
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Not only is there a significant discrepancy between our null result and the intensity 
calculated from the model, there is a factor of 20 difference between the observed O VI 
column density (~ 1.6 x 10^^ cm~^, Shelton & Cox (1994), Oegerle et al. (2002), see below) 
and the column density calculated from the model (2.7 x 10^^ cm~^), and a significant 
discrepancy between the 2 sigma upper limit on the observed C III 977 A intensity (7300 
photons cm^^ s^^ sr^^) and the published C III prediction (> 8500 photons cm^^ s^^ sr^\ 
Figure 2 Breitschwerdt (2001)). The number and magnitude of these discrepancies practically 
ehminate this class of models. 

Generalized Multi-Phase Model of the Local Bubble 

The lesson learned from evaluating the previous model was that a bubble consisting 
of a single, moderate temperature, overionized plasma is disallowed because it contains too 
many O VI ions and emits too many photons from its O VI and C III ions. Furthermore, a 
bubble consisting of a single, collisional ionizational equilibrium plasma which is hot enough 
to produce the observed soft X-ray surface brightness is disallowed because it contains too 
few O VI ions (see below). In contrast, suppose the Local Bubble contains multiple distinct 
plasmas, a collisional ionizational equilibrium, million degree. X-ray emissive plasma and 
a several hundred thousand degree, O Vl-rich plasma. Cool clouds are also known to he 
within the Local Bubble so such a topology is reasonable. One example of an allowed form 
is a hot bubble, bounded by transition temperature gas and possibly containing additional 
transition zones around embedded cool clouds. Here, we consider such a generalized multi- 
phase model, in order to determine if it can meet the observational constraints or if more 
unusual circumstances are required. 

The hotter component in this multi-phase model accounts for the observed soft X-rays, 
but few O VI photons. Based on the temperature found from a soft X-ray analysis (10^'^ K, as- 
suming collisional ionizational equilibrium, Kuntz & Snowden (2000)), the observed ROSAT 
R12 band countrate seen toward this part of the filament (590 xlO^® counts s^^ arcmin^^, K. 
D. Kuntz, private communication), the conversion between the ROSAT R12 band countrate 
and the emission measure (1 x 10~^ counts s"^ arcmin"^ = 7.15 x 10"^ cm~^ pc~^, Snowden 
et al. (1997)), the fraction of oxygen atoms in the O VI ionization level at this temperature 
(3.2 X 10~^, Mazzotta et al. (1998)), the solar abundance of oxygen (8.5 x 10~^, Grevesse & 
Anders (1989)), the ratio of free electrons to hydrogen nuclei (1.2), and the electron-impact 
excitation rate coefficient of ShuU & Slavin (1994), we estimate that the ~ 10^ K region 
produces only ~60 photons cm~^ s~^ sr~^ in the O VI resonance line doublet from a column 
density of only 5 x 10^^ O VI ions cm~^. 

The primary source of O VI ions and resonance line photons would be the transition 
temperature gas. Althought it need not be in collisional ionizational equilibrium, we will 
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assume its temperature is between 1 x 10^ K and 1 x 10^ K for the following calculations. 
The Local Bubble's column density of O VI ions is roughly 1.6 x 10^^ cm~^ (see following 
subsection). By assuming that the gas is in thermal pressure balance with the soft X-ray 
emitting plasma {Pth/k ~ 15,000 K cm^^, Snowden et al. (1998)), we estimate its electron 
density. The temperature is not well known, but most of the O VI ions are likely to be 
near their ionizational equilibrium temperature (3.2 x 10^ K), and the emission equation is 
relatively insensitive to temperature for 1 x 10^ ~ T ~ 1 x 10^ K. An O Vl-rich plasma, with 
the above column density, a temperature of 6.3 x 10^ K, and a thermal pressure of 15,000 K 
cm~^ emits 450 photons cm~^ s~^ sr~^, which is within the observational 2 sigma upper limit 
on the intensity. Similar plasmas with temperatures as low as ~ 1 x 10^ K emit even less 
intense O VI radiation. Thus, a two phase model could meet the observational constraints 
without requiring unobserved astrophysics or unusual conditions. 

Simulations of A Hot Bubble with Transition Zones 

In this section, we examine the multi-phase bubble simulations of Smith & Cox (2001) 
and the evaporating cloud simulations of Slavin (1989). In the models of Smith & Cox 
(2001), the Local Bubble resulted from two or three supernova explosions occurring up to 
several miUion years ago. The interior of the structure consists of hot (~ 10^ K), nearly 
collisional ionizational equilibrium plasma, while the periphery of the structure consists of 
an intermediate temperature transition zone. Although the simulations do not explicitly 
include the transition zones surrounding the embedded cool clouds, such zones arc thought 
to harbor O VI ions (Slavin 1989; Ocgcrlc ct al. 2002) and so should contribute to the flux 
of O VI resonance line photons. Here, we will use the simulations of Slavin (1989), which 
predict their intensity. 

After performing several hydrodynamical simulations. Smith & Cox (2001) decided that 
the best model lies intermediate between their various simulations. Their simulations yield 
predicted O VI intensities ranging from 190 to more than 9,500 photons cm~^ s~^ sr~^, (after 
the extra factor of 47r in their Figure 19 (Smith, private communication) is divided out). 
To this intensity must be added the intensity emitted by the transition zone on the local 
cloud and, possibly, transition zones on other cool clouds embedded in the Local Bubble. 
For the transition zone on the cool cloud surrounding the Sun, Slavin (1989) predicted an 
O VI doublet intensity of 250 photons cm~^ s~^ sr~^. The sum of the Smith & Cox (2001) 
and Slavin (1989) predictions (440 to > 9800 photons cm^^ s"^ sr~^) marginally overlaps 
the 2 sigma upper limit reported in this paper. If other cool clouds reside along the line of 
sight, then the predicted O VI intensity would be even higher. Thus, this particular model is 
heavily constrained by the null results. We would have lost confidence in this type of model 
had it not been possible to simultaneously satisfy the O VI emission, O VI column density. 
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and ROSAT sott X-ray constraints with the above generahzed multiphase model. 



5.2. Limits on Electron Density, Pressure, Pathlength, Cooling Timescale 

Here, we will combine the intensity upper limit with column density estimates taken 
from the literature in order to calculate the upper limits on the electron density and thermal 
pressure in the O Vl-rich plasma. The first data-derived estimate of the typical O VI column 
density of the Local Bubble (Nqvi ~ 1-6 x 10^^ cm^^ (Shelton & Cox 1994)) was found via 
a statistical analysis of several dozen Copernicus lines of sight, most of which terminated far 
beyond the Local Bubble. Later, FUSE observed a couple dozen lines of sight entirely within 
the Local Bubble. The work of Welsh et al (2002) and Oegerle et al. (2002) suggest that 
the O VI column density within the Local Bubble may be as little as a few xlO^^ cm~^ in 
some directions, but that the column density increases roughly with distance from the Sun 
and the total column density along a 100 pc radial path is consistent with values of 7 x 10^^ 
cm~^ to 1.6 X 10^^ cm~^. Here, we take the column density from the published survey (i.e. 
1.6 X 10^^ cm-2). 

Using the emission equation of ShuU & Slavin (1994), the above column density, the 
2 sigma upper limit on the doublet intensity, and the temperature at which O VI is most 
prevalent in coUisional ionizational equihbrium plasmas {T — 3.2 x 10^ K, Mazzotta et al. 
(1998)) yields an estimated upper hmit for the electron density, rie, of 0.021 cm~^. Thus, the 
estimated upper limit for the thermal pressure, from Pth/k — 1.92neT, is 13, 000 K cm~^. If 
the temperature were significantly higher than the assumed temperature, then the pressure 
would be greater. If it is not, then the total (thermal and nonthermal) pressure in the 
O Vl-rich region of the Local Bubble (presumably the transition zones) could be brought 
into balance with that of the X-ray emitting region of the Local Bubble (presumably the 
interior, Pth/k ~ 15,000 K cm~^) by assuming that the cooler, denser O Vl-rich regions have 
a greater non-thermal pressure than the hotter, more rarefied X-ray emitting region. 

Furthermore, the O Vl-rich material occupies a pathlength (A/) of at least ~5900A'"ov'//''^e 
in cm, (incorrectly type-set in Shelton (2002)), or ~ 0.15 pc, consistent with transition zones. 
The cooling timescale solely due to cooling through the O VI doublet is 

47r(sr)/oy7 Uovi 

where k is Boltzman's constant, riovii noxy, nn, and tit are the densities of O VI ions, oxygen 
atoms, hydrogen nuclei, and particles, respectively, sr means steradians, and Uovi/noxy — 0.24 
at 3.2 X 10^ K (Mazzotta et al. 1998). Thus, the coohng timescale due solely to O VI resonance 
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line emission is at least 2.0 x 10^ years. These Local Bubble values are tabulated in Table 4, 
along with the Galactic halo values calculated in the following subsection. 

5.3. Galactic Halo 

Previous FUSE studies of emission from O VI in the interstellar medium (ISM) outside of 
superbubbles or supernova remnants observed far longer paths (Shelton et al. (2001); Dixon, 
et al. (2001); Shelton (2002); and Welsh et al (2002)). In each of those cases, the instrument 

had been directed toward a high Galactic latitude direction {\b\ > 40°). Given that the 
local region contributes little O VI resonance line intensity, the vast majority of the intensity 
observed on the relatively unobscured, high latitude observations must originate beyond 
the Local Bubble. Extragalactic sources are ruled out by the velocities of the observed 
emission features and so this emission presumably resides in the Galactic thick disk and 
halo. Subtracting the Local Bubble's doublet intensity from the averages of the intensities 
observed by Shelton et al. (2001); Dixon, et al. (2001); Shelton (2002); and Welsh et al (2002) 
leaves 3500 to 4300 photons cm~^ s~^ sr~^, attributable to the Galactic halo. 

The column density through the halo is roughly 3.0 x 10^^ cm~^ (Savage et al. 2000). 
Thus, using the method of Section 5.2, we find the electron density in the halo's O Vl-rich gas 
to be 0.0050 to 0.0061 cm'^ The thermal pressure is 3100 to 3800 K cm"^, consistent with 
the range of total pressures in models of the Galactic thick disk and lower halo (Boulares 
& Cox 1990; Ferriere 1995). The pathlength is at least 110 pc. The cooling timescale is 7 to 
8 X 10^ years. 

6. Summary 

• The Local Bubble's O VI intensity has been examined through a shadowing observation 
directed toward / = 278.6°, b = —45.3°. No emission was seen. The tightest 2 sigma 
upper limits on the 1032 and 1038 A emission lines are 530 and 500 photons cm~^ 
g-i g2~-i^ respectively. Given the theoretical intensity ratio (2.0 to 1.0), the tightest 2 
sigma upper limit on the resonance line doublet is 800 photons cm~^ s~^ sr~^. 

• The class of models in which the Local Bubble is assumed to be a single, moderate 
temperature (~ 4 x 10'* K), extremely overionized plasma is virtually ruled out. The 2 
sigma upper limits on the O VI and C III intensities exceed the model predictions, as 
does the observed O VI column density. The class of models in which the Local Bubble 
is a single, hot (10^ K), coUisional ionizational equilibrium plasma are also ruled out. 
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Such a model could not produce the observed O VI column density. However, the 
combined O VI and ROSAT soft X-ray observations allow models in which the Local 
Bubble contains multiple, separate phases. The X-ray emission can be attributed to the 
hot zone and the O VI ions can be attributed to transition zones. Existing simulations 
of a multiphase bubble combined with simulations of a single evaporating cloud are 
tightly constrained by our 2 a upper limit on the doublet intensity. Considering that 
multiple clouds are expected to reside along the typical line of sight, the expected O VI 
intensity falls well below the upper limit, suggesting that transition zones differ from 
current models. 

• The estimated 2 sigma upper limit on the electron density in the Local Bubble's O VI- 
bearing plasma is 0.021 cm"^. If the temperature is near the coUisional ionizational 
equilibrium temperature for O VI (3.2 x 10^ K), then the thermal pressure is < 13, 000 K 
cm~^, which is weaker than that of the soft X-ray emissive portion of the Local Bubble. 
The discrepancy could be resolved if the O Vl-bearing gas were to have a greater than 
assumed temperature or a significant nonthermal pressure. 

• The timescale for cooling soley due to the O VI resonance line emission is proportional 
to the thermal energy in the gas and inversely proportional to the O VI intensity. For 
the Local Bubble, this timescale is > 2 x 10^ years. 

• The Galactic halo's O VI doublet intensity has been found by subtracting the Local 
Bubble result from other observations. It is 3500 to 4300 photons cm~^ s"^ sr~^. 

• The Galactic halo's intensity implies an electron density of ~ 0.0050 to ~ 0.0061 cm~^. 
If the temperature is near the coUisional ionizational equilibrium temperature for O VI 
then the thermal pressure is ~ 3100 to ~ 3800 K cm~^. This pressure is consistent 
with the total pressure estimates in models of the thick disk and lower halo. 

• The O VI in the Galactic halo cools on a timescale {t oc thermal energy / O VI intensity) 
of 7 to 8 X 10^ years. 

We thank the referee, Steve Snowden, for volunteering the DIRAS corrected IRAS and 
ROSAT images and critiquing the manuscript, K. D. Kuntz for providing the ROSAT mea- 
surements and commenting on the manuscript, B.-G. Andcrsson for discussing the telescope's 
specifications, and Tim Heckman for reviewing the manuscript. We acknowledge Welsh et 
al (2002), which compared observed and modeled C III intensities for other directions. This 
research has been supported by NASA grant NAG5-10394 (through the FUSE Guest Investi- 
gator program) and NASA grant NAG 5- 10807 (through the Long Term Space Astrophysics 
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Table 1. 



Program Id 


Start Date 


LiF lA Exposure Time 


Night Time 


LWRS Direction 


Number 




(UT) 


(ksec) 


(ksec) 


(l,b) 


12050501 


Aug 


■ 21, 


1999 


84.5 


27.8 


278.58°, - 


■45.31° 


12050510 


Aug 


• 24, 


1999 


17.7 


4.5 






12050601 


Aug 


. 27, 


1999 


18.2 


5.6 


278.59°, - 


45.30° 


B1290101 


Aug 


. 13, 


2001 


52.5 


9.6 


278.63°, - 


■45.31° 


B1290102 


Aug 


. 16, 


2001 


55.0 


15.4 






Total 








227.7 


62.8 


278.6°, - 


■45.3° 



Table 2. Upper Limits on the O VI Emission Lines 



Emission Region 


Intensity (ph s ^ cm ^ sr ^) 
Night Only 


Intensity (ph s ^ cm ^ 
Day+Night 


Combined Dataset 


1031.93 A 
1037.62 A 


180 ± 330 ± 14% 
-180 ± 310 ± 14% 


20 ± 230 ± 14% 
80 ± 210 lb 14% 


Tightest 1 Sigma Limits, inch 14% Systematic Uncertainty 


1031.93 A 

1037.62 A 


150 


280 


Tightest 2 sigma Limits, inch 14% Systematic Uncertainty 


1031.93 A 
1037.62 A 


500 


530 
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Tablc 3. Other Transitions 





I-? dct" \A/ Qiri^ 1 f^T\ cfi' h 
JTLCol V V dVclcii^ Lii 


liilCiiollj 


Notes 






{ Tr\\\r\'i'r\T\G /^m ^ c ^ ct* ^\ 






C I 


945.6 A 


-300 ± 390 


B 




C I 


1122.3 A 


-240 ± 370 


N 




C II 


1037.0 A 


30 ± 290 


N 




C III 


977.0 A 


4700 ± 1300 


N, 


SiC IB 






2600 ± 1000 


N, 


SiC 2A 


N I 


954.0 A 


-200 ± 740 


N 




N II 


916.7 A 


2300 ± 1100 


N, 


SiC IB 


N III 


991.6 A 


-170 ±850 


N 




Mg II 


946.7 A 


-430 ± 680 


N 




Si II 


1023.7 A 


-200 ± 190 


B 




S II 


906.9 A 


-240 ± 1100 


N 




S III 


1015.5 A 


170 ± 240 


B 




S IV 


1062.7 A 


-390 ± 300 


N 




S IV 


1073.0 A 


300 ± 280 


N 




S VI 


933.4 A 


170 ± 690 


N 




S VI 


944.5 A 


1000 ± 840 


N 




Fe II 


1144.9 A 


270 ± 220 


B 




Fc III 


1122.5 A 


-700 ± 350 


N 





The reported error bars reflect the 1 sigma random uncertainties. All measurements are 
subject to the net systematic uncertainty of 14%. "B" denotes both the day and night 
portions of the data and "N" denotes the satellite-night portion of the data. In cases where 
airglow emission or scattered light affected the satellite-day data for the spectral region of 
interest or the satellite-day spectrum undulated substantially, we used only the 
satcllitc-night portion of the data. We also used the night spectrum for the Mg II 
observation because it led to a more constraining upper Umit. For the measurements in this 
table, the tighter pulseheight cuts were used in processing the B12901 portion of the data. 
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Table 4. 



Region 


Tie 


Pth/k 


Al 


t 




(cm~^) 


(K cm-3) 


(pc) 




Local Bubble 


< 0.021 


< 1.3 X 10^ 


~ 0.15 


> 2 X 10^ 


Galactic Halo 


0.0050 to 0.0061 


3100 to 3800 


~ 110 


7 to 8 X 10^ 
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Fig. 1. — The Local Bubble (LB) is sketched in red in the center of the cartoon. The 
LB's size and shape were drawn from Snowden et al. (1998). The gray spots within the 
LB represent the cool clouds. The dotted, wiggly lines beyond the Local Bubble depict the 
outline of the Local Cavity, whose shape and size has been taken from the 50 mA Na I D2 
projected contours in Sfeir et al. (1999). To the right and lower right of the LB are cartoons 
representing the nearby Loop I superubble and the observed filament. Note that the filament 
is well beyond the calculated extent of the Local Bubble. An axis has been superimposed. 
The Sun lies at the origin. The Galactic center (GC) is to the right, the Galactic anticenter 
(GAG) is to the left, the north Galactic pole (NGP) is at the top, and the South Galactic 
pole (SGP) is at the bottom. Each of the tick marks is 200 pc from the center of the grid. 
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Fig. 2. — Right Panel: Map of the DIRBE-corrected IRAS 100 fim data for a portion of 
the sky centered on / ~ 280°, b ~ —45°. The "shadowing" filament discussed in this paper 
runs diagonally across the image. Left Panel: Map of the ROSAT I keV data for the same 
projection. Regions of missing data are black. By comparing these two maps, we can see that 
the "shadowing" filament corresponds with the sharply delineated arch of low ^ keV surface 
brightness in the soft X-ray map where the filament blocks soft X-ray emission from the 
Galactic halo. The finite "zero-point" surface brightness (~ 400 x 10~^ counts s~^ arcmin"^) 
is due to the Local Bubble emission. Both panels are adaptations of figures in Snowden et 
al. (1997). 
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Fig. 3. — Infrared 100/xm image of an opaque and bulbous section of the filament. The 
overlayed boxes indicate the placement of the FUSE LWRS aperture. The B1290101 and 
B1290102 observations point in the same direction but with slightly different angles. Their 
aperture location is labeled "A" and "B" on this image. At their location, the lOO/xm inten- 
sity is 7.2 MJy sr^^. The 12050601 pointing direction is labeled "C". For this observation, 
the location of the FUSE aperture crosses a pixel boundary in the lOO/zm map, such that 
the average 100/xm intensity is 7.2 MJy sr~^. The 12050510 and 12050501 pointing direc- 
tions are also nearly identical. They are labeled "D" and "E" on the map and have lOO/zm 
values of 7.0 MJy sr~^. The lowest lOO/xm intensity on this image is 2.0 MJy sr-^ Note 
that the filament's infrared brightness in excess of the background level implies a degree of 
obscuration that is in very good agreement with Penprase et al. (1998) 's quoted E(B-V). 
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B12901 

Night 



<^ 6x10 




VI VI 



Both 



^^^^^ 




6x10 



3x10^ - 



1020 1025 1030 1035 1040 1045 
CALFUSE Wavelength (A) 




1020 1025 1030 1035 1040 1045 
CALFUSE Wavelength (A) 



3x10 



6x10 



3x10 




1020 1025 1030 1035 1040 1045 
CALFUSE Wavelength (A) 



Fig. 4.— The 1020 to 1045 A region of the LiF lA spectra extracted from the B 12901 
(top panel), 120505 (middle panel) and 120506 (bottom panel) data. Each panel presents 
the spectra drawn from the sateUite- night portion of the data (raised by 6 x 10^ pho- 
tons s~^ cm""^ sr~^ A~^), the spectrum drawn from the satellite-day portion of the data 



(raised by 3 x 10^ photons s ^ cm ^ sr' 



A~^), and the spectrum drawn from the combined 
day and night data (unadjusted). The wavelength scales have not yet been corrected. Each 
spectrum has been binned by 11 pixels The Earth's atmospheric airglow lines are marked 
with © symbols. The wavelengths of the O VI resonance lines and C II 3s^ S1/2 to 2p^ P3/2 
line are marked on the top and middle panels. No O VI or C II emission is evident. 



Night 



Doy 



-24- 



4x10''" 7 



Satellite Night 



E 
o 



2x10^- 



Q_ 



1 



% 1x10^^ 

-I— ' 

c 





•< 



E 
o 



4 - 

4x10^ r 



3x10 



4 - 



2x10 



Q_ 



I 1x10^7 
c 



Day + Night 






1028 



1038 1040 



1030 1032 1034 1036 
LSR Wavelength (A) 

Fig. 5. — Co-added spectrum from the night-only portion of the data (top panel) and 
co-added spectrum from both day and night portions of the data (bottom panel). The 
wavelength scale has been corrected and shifted to the LSR reference frame. The spectral 
data has been binned by 11 pixels. The night spectrum is flat, except for variations due to 
noise. No excess intensity around 1032, 1038, or 1037 A is observed. Neither is there a 1032 A 
feature in the "Day+Night" spectrum. The sharp thin emission feature near 1031.5 A is only 
present during the satellite-day portion of the data and so is not attributable to interstellar 
material. Similar features have appeared in previous daytime blank sky observations (Shelton 
et al. 2001). The "Day-I- Night" spectrum wavers significantly between 1036 and 1038 A, 
preventing us from searching for the 1037 and 1038 A hues. 



